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Abstract 



We develop a neutrino transport theory, in terms of the real-time non- 
equilibrium Green's functions, which is applicable to physical conditions ar- 
bitrary far from thermal equilibrium. We compute the coherent neutrino 
radiation in cores of supernovae by evaluating the two-particle-two-hole (2p- 
2h) polarization function with dressed propagators. The propagator dressing 
is carried out in the particle-particle channel to all orders in the interaction. 
We show that at two loops there are two distinct sources of coherence effects 
in the bremsstrahlung. One is the generically off-shell intermediate state 
propagation, which leads to the Landau-Pomeranchuk-Migdal type suppres- 
sion of radiation. We extend previous perturbative results, obtained in the 
leading order in quasiparticle width, by deriving the exact non-perturbative 
expression. A new contribution due to off-shell finial/initial baryon states is 
treated in the leading order in the quasiparticle width. The latter contribu- 
tion corresponds to processes of higher order than second order in the virial 
expansion in the number of quasiparticles. At 2p-2h level, the time compo- 
nent of the polarization tensor for the vector transitions vanishes identically 
in the soft neutrino approximation. Vector current thereby is conserved. The 
contraction of the neutral axial vector current with tensor interaction among 
the baryons leads to a non-vanishing contribution to the bremsstrahlung rate. 
These rates are evaluated numerically for finite temperature pure neutron 
matter at and above the nuclear saturation density. 
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I. INTRODUCTION 



Neutrino production in baryon encounters is among the fundamental processes by which 
compact stars lose their energy. The reactions can be arranged, in general, according to 
the number of participating baryons, as the phase space arguments play the central role in 
controlling their temperature and density dependence ■ In the case of neutrino pair 
bremsstrahlung, the leading order process in the density virial expansion is the two-body 
reaction 



B x + B 2 -»• B x + B 2 + u f + v f , (1) 

where B stands for a baryon, vj (pf) for a neutrino (anti-neutrino) of flavor / = e, /i, r. Note 
that the subleading order process (i.e. the one in the absence of the spectator) vanishes for 
identical particles, as an on-shell propagating particle cannot radiate. 

The matter in neutron stars is highly degenerate for temperatures typically below a 
MeV and the elementary excitations are quasiparticles with well-defined energy-momentum 
relation. Produced neutrinos are typically "soft" with energies of order of temperature. 
In this limit the intermediate quasiparticle propagator diverges as 1/u and the amplitudes 
of the neutrino absorption, scattering, and emission turn out formally divergent as l/u> 2 . 
The infrared behaviour of the in-medium rates, however, is dominated by the neutrino 
phase space factors, rather than the infrared divergence of the amplitudes and the rates of 
the bremsstrahlung and its space-like analogues remain finite. At the same time, at low 
temperatures, the contribution from the infrared region to the rate of the bremsstrahlung 
is negligible. The combined effect of the cancellation of the infrared divergence and the 
vanishing contribution from the low frequency region makes the quasiparticle approximation 
to eq. (0) applicable in cold neutron stars. 

During the first several tens of seconds after a supernova explosion and core collapse 
the temperature of the dense nuclear matter is of the order of several tens of MeV. The 
neutrino bremsstrahlung is then suppressed, because the formation length of the neutrino 
radiation is of the same order of magnitude as the mean free path of a baryon p]-|TT| . The 



collective effects become important on the radiation scale (i.e. the role of the spectator in the 
reaction ([I]) is taken over by the medium) because the baryon undergoes multiple scattering 
during the radiation. The underlying mechanism is the Landau- Pomeranchuk-Migdal (LPM) 
quenching of the radiation, first introduced in the context of QED |IT|. The central role 



in the theory is played by the formation length of radiation If. If the mean-free-path of a 
baryon is much larger that the formation length Z m f p ^> If then the radiation reduces to 
a sum of separate radiation events, each of which is well described by the Bethe-Heitler 
spectrum. In the opposite limit Z m f p <C If the individual scattering events are unresolved 
and the radiation spectrum takes the Bethe-Heitler form for a single scattering event. In 
the intermediate regime, when l m f p ~ If, the radiation amplitudes for scattering off various 
centers interfere destructively and the radiation is suppressed (Landau- Pomeranchuk-Migdal 
effect; for a review see refs. P,|12| ) . 

The rates of neutrino-nucleon processes are commonly expressed through phase space 
integrals over the contraction of the weak currents with the polarization function of the 
nuclear medium. The polarization function (or structure function) of the supernova/neutron 
star matter has been subject of many studies P,|T3|-|TE|] . The modifications of reactions rates 
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by the spatial correlations among (on-shell) quasiparticles have been studied within the 



Fermi-liquid theory [fL3|l , the one-boson exchange interaction theory the relativistic 



random phase approximation the variational approach [16[], and combinations thereof 



l7| . [T8|| . The spatial correlations tend to suppress reaction rates in general, although their 



impact on the supernova physics is model dependent fl3|-|i"8|. 

The common strategy of incorporating the LPM effect in the neutrino-nucleon interaction 
processes is to add a quasiparticle damping in the intermediate state propagator by replacing 
uj by uj + ry |7|-|| . In the soft neutrino limit the vector current coupling does not contribute 
by virtue of the vector current conservation (CVC) and the net contribution comes from 
the axial- vector transitions via baryon spin-flip. The above modification of the intermediate 
state propagator then leads to an ansatz for the nucleon spin structure function: S a oc 
7o-/(u; 2 + 7^) [0-0, where 7 CT is the nucleon spin-flip collision rate. The ansatz generalizes 
the quasiparticle picture, in a semi-phenomenological manner, by including the temporal 
correlations among the quasiparticles in the leading order in the quasiparticle width. The 
microscopic justification of this phenomenology emerges from the various formulations of the 
finite temperature quantum filed theory, e. g. the thermo-field dynamics [|TJ]] or the closed 
diagram formalism in the Schwinger-Keldysh technique [0. A microscopic computation is 
not straightforward, however. For example, the polarization function of the medium can be 
computed at one-loop, including the quasiparticle width to all orders in 7, however a priori 
the current conservation is not guaranteed at this level. The reason, in part, is that the 
"more complicated" higher order in loop expansion diagrams contribute at the same order 
as the single loop Ul2| . 

In a previous paper we carried out a microscopic computation of the bremsstrahlung, 
including the LPM effect, at the one-loop level in a formalism based on the quasiclassical 
Kadanoff-Baym transport equation [^DJ. Here we extend this computation to two- loops and 
partially modify our approach to include the propagator and vertex renormalization on the 
same footing and to including the tensor force explicitly. The extension to two-loops is 
motivated by the following. The long range phenomena, driven by the weaker attractive 
part of the baryon-baryon interaction, are sensitive to the resummation in the particle-hole 
{ph) channel. On the other hand, as well known, one should fully resum the particle- 
particle (pp) channel to treat the hard core of the baryon-baryon interaction. Therefore, 
the ph channel can be treated perturbatively by expanding in the number of particle-hole 
loops, while the pp channel must be treated non-perturbatively by a full resummation of 
the ladder diagrams. Thus, the separation of the long range and short range phenomena 
dictates the manner in which the diagrammatic expansion is carried out. The dressing of 
the single particle propagators occurs in both channels and can be treated either explicitly, 
say, by considering higher order self-energies attached to a propagator, or, alternatively, by 
condensing it in the width of the propagator spectral function. As a consequence of the 
separation of the scales, the short-range correlations can be condensed in the propagator 
width on the scales relevant for the long-range phenomena. The imaginary part of a single 
loop in the ph channel vanishes in the time-like region of the phase space, which is relevant for 
the particle production. A finite result emerges when one dresses the propagators by either 
extending the resummation in the ph channel to two and higher loops and/or by dressing 
the propagators in the pp channel to all orders. Ignoring the latter resummation, i.e. using 
the quasiparticle propagators in the two-loop expansion, misses a number of short-range 
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collective effects, such as the LPM quenching of the radiation due to multiple scattering. 
On the other hand, summing only the ladders in the pp channel does not recover the vector 
current conservation in the radiation process (in, at least, a transparent manner). Therefore 
a natural choice, motivated by the separation of the short and long range phenomena, is 
to truncate the ph channel at two-loops and to resum the pp channel to all orders. The 
situation is reminiscent of the parquet resummation scheme in the first iteration, where in 
both channels the driving force is the bare baryon-baryon interaction. 

Early studies of the bremsstrahlung at the quasiparticle level modelled the strong force 
using the T-matrix interaction ||, the free-space one-boson exchange interaction || and 
their in-medium modifications H supplemented with a hard core modelled in the spirit 
of the Fermi-liquid theory. The explicit use of the tensor interaction turned out to be 
crucial as there are significant cancellations among different diagrams, and the surviving 
contribution is due to a non-trivial contraction between the operator structures of the weak 
and strong interactions (tensor force) [|[]. This motivates our ansatz for the driving force 
in the particle-hole (ph) channel of nuclear interaction, which includes explicitly the tensor 
force contribution. We do not attempt, in the present work, to go beyond the one-pion 
exchange approximation for several reasons, one being that the non-perturbative treatment 
of the interaction does not change the spin, isospin, and tensor operator structure of the 
interaction, and important cancellations in the radiation matrix elements will be preserved 
in a more advanced treatment. We also want to be able to isolate the finite width effects in 
our comparisons to the earlier work done in the one-pion exchange approximation [f|||. The 
situation is different in the particle-particle (pp) channel, where the short-range correlations 
have to be treated in a non-perturbative manner by summing up the ladder diagrams to all 
orders. We do this in the finite-temperature Brueckner theory. 

The paper is organized as follows. In Section 2, starting from the Kadanoff-Baym formal- 
ism, we derive a single-time transport equation for (anti)-neutrinos with collision integrals 
driven by (anti)-neutrino coupling to baryons via the polarization tensor of the medium. 
The polarization tensor is computed in the 2p-2h approximation in Section 3. The sum- 
mation of the ladder diagrams in the pp channel within the finite temperature Brueckner 
theory is described in Section 4. Section 5 evaluates the phase space integrals and neutrino 
bremsstrahlung emissivities. The numerical results are presented in Section 6. Section 7 
summarizes our main results. 



II. NEUTRINO TRANSPORT FORMALISM 

A. Neutrino propagators 

The theory of neutrino radiation can be conveniently formulated in terms of the real- 
time quantum neutrino transport. Let us start by defining the various time-ordered Greens 
functions of massless Dirac neutrinos. These can be written in the generic matrix form 

• q _ .( s c 12 s< \ _ ( (TiKztWfaj) - (feMa*)) \ _.( s u - S u + \ 

l ^ - 1 \ s> s« 12 ) - \ (^(xMx 2 )) (m Xl mx 2 )) stf s+ + ) ' {I) 

where tp(x) are the neutrino field operators, ip = 7 C V% T is the chronological time ordering 
operator, and T is the anti-chronological time ordering operator; the indexes 1 = Xi, 2 — 
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x 2y .. collectively denote the space-time and discrete quantum numbers. The neutrino matrix 
propagator is further assumed to obey the Dyson equation, 



Sfa, x 2 ) = Sj (x 1 ,x 2 ) + Sofa, x 3 )fl{x 3 , x 2 )S{x 2 , x x ) 

= Sofa, x 2 ) + Sfa, x 3 )tt(x 3 , x 2 )S (x 2 , xi), (3) 

where Sofa,x 2 ) is the free neutrino propagator and Sq (xi, x 2 )Sofa, x 2 ) = a z 5{x\ — x 2 ), 
a z is the third component of the Pauli matrix, Q is the neutrino proper self-energy and 
we assume integration (summation) over the repeated variables. The self-energy Q is a 
2x2 matrix with elements defined on the contour in terms of the Dyson equation. The 
quasiclassical neutrino transport equation follows from the Dyson equation in the 'conjugate 



subtracted' form 21,22 : 



iSfa,x 2 ) $ X2 -i $ Xl Sfa,x 2 ) = Sfa,x 3 )Qfa,x 2 )q^-o JL Qfa,X3)Sfa,x 2 ), (4) 

Note that the initial correlations are neglected in eq. (|j). The set of the four Green's 
functions above can be supplemented by the retarded and advanced Green's functions which 
are defined as 

is* = eih-h^Ufai^fa)}), is* = -e(h-h)(Ufa)rffa)\), (5) 



where 8(x) is the Heaviside step function on the real-time contour defined as d9(x)/dx = 
a z 8(x). The retarded and advanced Green's functions obey integral equations in the quasi- 
classical limit. The relations between the six Green's functions are listed in the Appendix A. 
The transport equation for the off-diagonal elements of the matrix Green's function reads 



d X3 -Ke^(x 1 ,x 3 ),^ >1< (x 3 ,x 2 ) - ^eS R {x 1 ,x 3 ),Q > ' < (x 3 ,x 2 ) 



= i {S>'<fa, x 3 ), n>><fa, %2)} + l {^' K fa, x 3 ),S>><(x 3 , x 2 )} , (6) 

where [ , ] and { , } stand for commutator and anti-commutator, respectively. In arriving at 
eq. (H) w e assumed the existence of the Lehmann representation for the neutrino propaga- 
tors; as a results we have 9te S R = Ke S A = Ke S and Ke £l R = Ke £l A = Ke SI. 

For the present purposes the neutrino dynamics can be treated semiclassically, by sep- 
arating the slowly varying center-of-mass coordinates from the rapidly varying relative co- 
ordinates. Carrying out a Fourier transform with respect to the relative coordinates and 
keeping the first-order gradients in the slow variable we arrive at a quasiclassical neutrino 
transport equation 

% {MeS-\q, x), S^iq, x)} pj? + i {3?e S(q, x), n>><(q, x)} pB . 

= S > ' < {q, x)Q > ' < {q, x) + fi >5< (g, x)S > ' < (q, x), (7) 

where q = (q, go) and x are the neutrino four momentum and the center-of-mass space-time 
coordinate, respectively, {. . .}p.b. is the four-dimensional Poisson bracket. The l.h.s. of eq. 
(0) is the precursor of the drift term of the Boltzmann equation. The second Poisson bracket, 
however, does not fit in the Boltzmann description and can be eliminated by an expansion 
of the neutrino propagator in the leading (quasi-particle) and next-to-leading order terms 
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in the small neutrino damping: S > ' < (q, x) = Sq' (q,x) + S x ' (q, x). A direct evaluation of 
the Poisson brackets decouples the l.h.s. of transport equation (|7|) to the leading order with 
respect to the small damping of neutrino/anti-neutrino states (Qrnf2(g, x)/dteQ(q, x) <C 1). 
The quasiparticle part of the transport equation 

i {®eS-\q, x), S^q, x)} pB = S>'<{q, x)V>><{q, x) + fi>><(g, x)S > ' < {q, x) (8) 

describes the evolution of the distribution function (Wigner function) of on-shell excitations 
with the l.h.s. corresponding to the drift term of the Boltzmann equation. The r.h.s. 
corresponds to the collision integral with the self-energies Q >,< (q,x) having the meaning of 
the collision rates. The advantage of this form of the (generalized) collision integral is that 
it admits systematic approximations in terms of the Feynman perturbation theory. The 
remainder part of the transport equation 

i {MeS-\q, x), S>'<(q, x)} pj? + i {Ke S(q, x), fi>'<(g, x)} PB = 0, (9) 

relates the finite width part of the neutrino propagator to the self-energies in a form of a 
local functional which depends on the local (anti-)neutrino particle distribution function and 
their coupling to the matter. 



B. On-shell neutrino approximation 



The on-mass-shell neutrino propagator is related to the single-time distribution functions 
(Wigner functions) of neutrinos and anti-neutrinos, f u (q,x) and f„(q,x), via the ansatz 



So(q,x) 



ITX i 



S (g - w u (q)) f u (q, x) - S (q + u u (q)) (1 - fp(-q, x)) 



(10) 



where oo u (q) = c\q\ is the on-mass-shell neutrino/anti- neutrino energy. Note that the ansatz 
includes simultaneously the neutrino particle states and anti-neutrino hole states, which 
propagate in, say, positive time direction. Similarly, the on-shell propagator 



So(q,x) 



VK 



u v {q) 



S (g - Mu(q)) (1 - fu(q, x)) - 5 (g + u v (q)) f„(-q, x) 



corresponds to the states propagating in the reversed time direction and, hence, includes 
the anti-neutrino particle states and neutrino hole states. 

To recover the Boltzmann drift term in the on-shell limit, we take the trace on both 
sides of the transport equation (0) and integrate over the ( ant i-) neutrino energy q Q . The 
first term on l.h.s. of eq. (0) reduces then to the drift term of the Boltzmann equation. The 
single time Boltzmann equation (hereafter BE) for neutrinos is obtained after integrating 
over the positive energy range: 

d t + d q u u (q)d x ] f u (q,x) = r ^Tr [0%, x)S>(g, x) - 0>(q, x)S<(q, x)] ; (12) 
J JO Ztc 

a similar equation follows for the anti-neutrinos if one integrates in eq. (0) over the range 
[-oo,0]. 
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Fig 1: The neutrino Dyson equation in terms of the Feynman diagrams. The dashed curve corresponds to 
the S'-propagator, which includes the neutrinos and anti-neutrino holes moving in the same time direction; 
(reverting the time-direction one finds the Dyson equation for anti- neutrinos and neutrino holes). The 
shaded loop is the baryon polarization tensor. The wavy lines correspond to the W ± , Z° boson propagators. 

The different energy integration limits select from the r.h.s. of the transport equations the 
processes leading to modifications of the distribution functions of (anti-)neutrinos. The sep- 
aration of the transport equation into neutrino and anti-neutrino parts is arbitrary, however 
is motivated by the observation that the fundamental quantities of neutrino radiative trans- 
port, as the energy densities or neutrino fluxes, can be obtained by taking the appropriate 
moments of BEs. These quantities are not symmetric with respect to the neutrino/anti- 
neutrino populations in general. E.g. the neutrino emissivities (energy output per unit time 
per unit volume) for processes based on /9-decay reactions are given by the zeroth order mo- 
ment of the anti-neutrino BE, and it is sufficient to consider only the BE for anti-neutrinos. 
In the case of the bremsstrahlung we have to eventually sum these equations; still the rela- 
tion of the transport self-energies to particular processes becomes transparent if one treats 
the transport equations separately. 



C. Collision integrals 

We adopt the standard model for the description of the neutrino-baryon interactions and 
write the neutral current interaction Hamiltonian in the from: 

H int = ^=T H T L , T H = ^(c v -c Al5 )<f>, r L =^(l- 75 )^, (13) 

where G is the weak coupling constant, ip and <fi are the neutrino and baryon field opera- 
tors, cy and ca are the dimensionless weak neutral-current vector and axial vector coupling 
constants. 

The diagrams contributing to the neutrino emission rates can be arranged in a perturba- 
tion expansion with respect to the weak interaction. The lowest order in the weak interaction 
Feynman diagrams which contribute to scattering, emission, and absorption processes are 
shown in the Fig. 1. The corresponding transport self-energies are read-off from the diagram 

-in>'<( qi ,x) = J ^ f ? 4 (27r) 4 ^ 4 ( gi - ft - qm q iS<( q2 ,x)iTi\mli<(q,x), (14) 

where II f (q) are the off-diagonal elements of the matrix of the baryon polarization tensor, 
rl q is the weak interaction vertex. The contact interaction ( |T3|) can be used for the energy- 
momentum transfers much smaller than the vector boson mass, q <C m^m^y. Let us first 
concentrate on the BE for neutrinos. Define the loss and gain terms of the collision integral 

as: 
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r>,<, 



q, x) 



00 rfgo Tr 
2tt r 



^(g,^^) 



(15) 



Substituting the self-energies and the propagators in the collision integrals we find for, e.g., 
the gain part: 



dq 



10 



2tt 



Tr< 



d 4 q d 4 q2 
(2vr) 4 (2vr) 4 



(2tt) - ft - g)rg- 



7T02 



M<?2j 



5 (g 2 - w v (g 2 )) fv(q2, 



$ (qo2 + w v (q 2 )) (1 - /s(-«2, &)) 



rt A ^ \ £ (fro - w v (qfi)) (1 - /„(fr , x)) n> A (g, 2) 



(If 



The loss term is obtained by replacing in eq. (|16|) the neutrino Wigner functions by the 
neutrino- hole functions f u (q, x) — > (1 — /„(g, a;)) and the anti-neutrino-hole Wigner functions 
by the anti-neutrino functions (1 — fo(—q, x)) — > fo(q, x). The terms proportional (1 — f v )f v 
and (1 — f v ){l— fp) in the gain part of the collision integral, I„(q), correspond to the neutrino 
scattering- in and emission contributions, respectively. The terms proportional f v (l — /„) 
and f u f p in the loss part of the collision integral, I„(q), are the neutrino scattering-out and 
absorption contributions. 

The loss and gain collision integrals for the anti-neutrinos can be defined in a manner, 
similar to the case of neutrinos, with the energy integration spanning the negative energy 
range 

r0 dq 0r 



r>,<, 



q,x) 



2tt 



-Tr 



n>'<(q,x)S^ < (q,x) 



(17) 



Using the above expressions for the self-energy and the propagators, we find, e.g., for the 
gain term: 



dq 



10 



2tt 



Tr< 



d 4 q d% fn _^x± f „ „n-pm 



00 (2tt) 4 (2tt)' 



-(2vr)V(g 1 -g 2 -g)r^ 



^(92) 



8 (g 2 - u u (q 2 )) f v {q 2 , x 



S (q 2 + u u (q 2 )) (1 - fp{-q 2 , x)) 



4A 



7T 



Uv{qi 



-S (q 10 + u v {qi)) fu(-qi,x)U> x (q, x) 



The loss term is obtained by making replacements in eq. (|18D analogous to those applied to 
eq. ([TBI) . The terms proportional f u f p and /p(l — fo) in the gain part of the collision integral, 
Ip(q), then correspond to the neutrino absorption and scattering-out contributions. The 
terms proportional (1 — /p)(l — f v ) and (1 — f&)f& in the loss part of the collision integral, 
Ip(q), are the neutrino emission and scattering-in contributions, respectively. Note that, 
when the neutrinos are in a thermal equilibrium with the baryons, the collision integrals for 
the scattering- in/scattering-out and for the absorption/emission cancel. Under the condi- 
tions of detailed balance the (anti-)neutrino distribution function reduces to the Fermi-Dirac 
form. 



D. Bremsstrahlung emissivity 

The neutrino-pair emissivity (the power of the energy radiated per volume unit) is ob- 
tained by multiplying the left-hand-sides of the neutrino and anti-neutrino by their energies, 
respectively, summing the BEs, and integrating over a phase space element: 



S 



jtlwf [Uq) + Mq)] ^ = Iw)~ 3 [I "' cm{q) ~ /r ' cm(<7)] (19) 



where in the collision integrals we kept only the terms which correspond to the processes 
with the neutrino and anti-neutrino in the final state (bremsstrahlung) 



/ 



d3( ll r>,<,em/ \ / \ f rf V <PQ2 (Pq l :i 



(2tt)3 v v " iy v ^> J (27r)32^(g 1 )(2 7 r)32^(g 2 )(27r) 4, 

5(^(90+^(92) - 9o )^(90 [1 - /»K(<?0)] [1 " /p(^(9 2 ))] A M (gi,g 2 )n>j v < (g,x), (20) 



and A mA = Tr [7^(1 — 7 5 ) ^7^(1 — 7 5 ) ^ 2 ]- The collision integrals for neutrinos and anti- 
neutrinos can be combined if one uses the identities n< A (g) = (— q) = 2ig B (g )SmII^(g); 
here gB^Qo) is the Bose distribution function and H^ x (q) is the retarded component of the 
polarization tensor. With these modifications the neutrino-pair bremsstrahlung emissivity 
becomes 

d 3 q 2 f d 3 qi f d A q 




(27r)32c^(9 2 ) J (2tt)32^(9i) J (2vr) 4 
(2tt) 4 <5 3 (9i + q 2 - q)8{u v {q l ) + u v (q 2 ) - q ) [w„(qi) + ^(92)] 

9B(qo) [1 - fuMqi))] [1 " /,K(q 2 ))] A M (gi,g 2 )SmnJ A (g). (21) 

We note that eq. ([H]) is applicable for arbitrary deviation from equilibrium, as the equilib- 
rium properties of the neutrinos and baryons have not been used in the derivation (e.g. the 
temperature of the bath drops out if one assumes an initially uncorrelated state). Therefore 
eq. fl2~T|) is applicable beyond the boundaries of the linear response theory or the S'-matrix 
theory which explicitly resort to the equilibrium properties of the system as a reference 
point. 



III. TWO-LOOP BARYON POLARIZATION FUNCTION 

In this section we start the implementation of the perturbative scheme motivated in the 
introduction. Our strategy is the separation of the long and short range phenomena in the 
ph and pp channels. Here we carry out the first step by expanding the particle-hole channel 
and truncating it at two loops. This fixes the amount of the long-range correlations in the 
theory. The short-range effects are condensed in the width of the particle-hole propagators, 
which is specified in a later section by summing the ladder diagrams. 



A. Baryon propagators 

Although we shall treat the baryon sector in the equilibrium limit, it is still useful to 
define the six Green's functions of the non-equilibrium theory, as in the case of neutrinos. 
The matrix Green's function of non-relativistic baryons is defined in the standard way 

iC i ( G ™ G< \ _ / (T0(xO0t(x 2 )) - (0t(x 2 )0(xO) \ _ / G- l2 G- l2 + \ 
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where 4>(x) are the baryon field operators. In terms of these operators the retarded and 
advanced function are defined as 



iG R 2 = e(t 1 -t 2 )({<P(x 1 ),^(x 2 )}), iG? 2 = -9{t 2 -t 1 )({^ 1 ),0 t (^)}). (23) 



The structure of the proper self-energy matrix E is identical to eq. Q2"2]) and its elements are 
defined via the Dyson equation for baryons: 

G(x 1 ,x 2 ) = Qo(xi, x 2 ) + Go(xi, x 3 )E(x 3 , x 2 )G(x 2 , x x ) 

= G Q (x 1 ,x 2 ) + G(xi,x 3 )E(x 3 , x 2 )G (x 2 , xi). (24) 

In a complete analogy to the neutrino sector, we approximate the Green's functions by their 
quasiclassical counterparts by defining center-of-mass and relative space-time coordinates 
and Fourier transform with respect to the relative space-time coordinates. In the equilibrium 
limit the dependence of the quasiclassical Green's functions on their center-of-mass space- 
time coordinate is trivial and can be dropped. The distribution function of the baryons is 
related to the off-diagonal elements of the matrix Green function by the exact relations 

- iG<(p) = a(p)f N (p), iG>(p) = a(p)[l - f N (p)\, (25) 

where a{p) = i[G R (p) — G A {p)} = i[G > (p) — G < (p)} is the baryon spectral function, /n(p) = 
[exp(/3(u — /i)) + l] _1 is the Fermi-Dirac distribution function, j3 = T _1 is the inverse temper- 
ature and fi is the chemical potential (relations (|25|) will be refereed to as the Kadanoff-Baym 
ansatz in the following). The quasiparticle energy, e p = p 2 /2m + ^teT l R (p)\ UJ=Ep follows from 

the solution of the Dyson equation G R (p) = uj — e p + i^smS R (p) . When damping of 

quasiparticle states is small, QmL R (p) 9fte£ R (p), the propagators can be decomposed 
into quasiparticle and background contributions, e.g., 

G<(p) ~ 2mz(p)f N (p)8(u - e p ) - E<(p) - + 0( 7 2 )- (26) 

[UJ E p ) 

Note that the self-energy appearing in the denominator of the second term of eq. fl2~6|) via 
the dispersion relation is restricted, to the leading order in damping, to the mass-shell. In 
equilibrium, 

iE<(p)=7(p)/ JV (p) J -iS > (p)= T (p)[l-/ JV (p)], (27) 

where "f(p) = — 29mE(p) is the width of the baryon spectral function. The wave-function 
renormalization, z(p), in the same approximation is 

^3mE(a/,p) - / , (28) 

where we used the integro-differential form of the Kramers-Kronig relation: 

= /— 3mS(a/,p) V (29) 
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On inserting the expression of the wave-function renormalization (|28|) in the expansion (|26] 
we find the final form of the propagator 



/dtu' "P 
T- ^T~, \2 - ^ - ^ - ^ M"). (30) 



da/ , „ P 

Note that this form of propagator renders the strict fulfillment of the spectral sum rule, 

/^a(p) = l, (31) 

at any order in the expansion with respect to the damping. 

Using the linear relations among the propagators, listed in Appendix A, we find for the 
causal propagator: 

G —^_ w - (e p + sfte£(» - p) 



[to - (e p + 3te£(p) - fx)} 2 + [3mE(p)] 2 



^ mS(p) r tanh ( ^] , (32) 



[u - (e p + 3te£(p) - + [9mS(p)]' 

where tanh (u>/2) = [1 — 2/jv(a?)] and e p = p 2 /2m. As the evaluation of the baryon polariza- 
tion function requires the causal and acausal Green's functions of the type G (q + p), we 
note here that, the denominator of such a function can be expanded in the limit vq <C u, 
where v <C 1 is the characteristic velocity of a baryon, 

d 

{u + e p ) - £p> + g- ~ oj - p ■ q/m - q— 3te£(p) - e q ~ u, (33) 

to the leading order. The approximation (|33|) will be referred in the following as the soft- 
neutrino approximation. We also employed the non-relativistic limit for baryons. If we use 
the ansatz j(— uS) = j(oo), which is exact in the phenomenological Fermi-liquid theory and 
will be verified in our microscopic calculations, then 

G-{±u,p) = ±- 2 =F t 2 7( "; P)/ L tanh ( ^] , (34) 

V yJ cu 2 + 7 (u;,p)V4 ^ u 2 + 1 {uo,p) 2 /A \2) 

-G++(±w,p) = ± - " ^- ± i 2 1{U ] P)/ 1 IA tanh ( ^) , (35) 

V ; uo 2 + 1 (oo,p) 2 /A u 2 + 7 (u;,p) 2 /4 V 2 / 

where the second equation follows from the relation [G (p)]* = —G ++ (p), valid in the 
momentum representation (see Appendix A). Thus both propagators are odd under the 
exchange of the sign of a;, a property which will be important in establishing the vector 
current conservation in the radiation processes discussed below. Since the dependence of 
the the quasiparticle width on the momentum is weak in the density and temperature range 
of interest it is useful to define momentum average quasiparticle width which a function only 
of the frequency. This approximation is implemented in the phase space integrations below. 



11 



B. The interactions 



The central ingredient of a bremsstrahlung process is the modelling of strong the inter- 
action. For the particle-hole interaction a reasonable, but not unique, choice is the one-pion 
exchange interaction combined with a contact interaction in the spirit of the Fermi-liquid 
theory: 

W*) = (—) ■ k ) D ~~( k ) • k) + f + h{tr ■ <r), (36) 

where f n is the pion decay constant, is the pion mass, D (fc) is the one-shell causal 
pion propagator, f and f\ are the coupling parameters of the Fermi-liquid theory, cr is 
the vector of the Pauli matrices. The non-relativistic reduction of the neutrino-neutron 
interaction vertex ( O ) is 

G 

r f = (^o - SA^^i) , (37) 

where Qa = 1-25 is the axial-vector coupling constant. 



C. Direct contribution to the polarization function 

The three topologically different direct diagrams (i.e. those which do not involve an 
exchange of outgoing particles) are shown in Fig. 2a-c. 







Fig. 2: The Feynman diagrams for neutrino- nucleoli interaction in the 2p-2h approximation. The vertical 
dashed lines correspond to the baryon-baryon interaction and the wavy lines to the Z° vector bosons. 
Exchange diagrams are shown below in Fig. 3. 
The analytical expression, corresponding to the Fig. 2a, is 



ill 



%) = /n 



(27T)' 



dk 



{27i) 8 S(q + p 4 - k - P3 )5(k + p 2 -pi)Tr \v(k)G~ + ( Pl )V(k)G + - (p 2 



Ti\r^G--{q + p 4 )V{k)D--{k)G- + {pz)V{k)D ++ {k)G ++ {q + p i )T u G + -{p 4 , 

where V(k) is the strong interaction vertex, which can be read-off from eq. (RBI) . The 
contribution of this diagram is readily recognized as a propagator dressing in the ph channel 
by means of a self-energy corresponding to an excitation of a single particle-hole collective 
mode. The analytical expression, corresponding to the Fig. 2b, is 



(3* 
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in 



-+,6 



(?) = n 



, 4 



i=l 



d 4 Pi 
(2tt)< 



dk 



{2n) 8 5(q + p A - k - p 3 )5(k + p 2 - Pl )Tr \v(k)G- + ( Pl )V(k)G + -(p 2 ) 



Ti[Y,G--{q + p A )V{k)D--{k)G- + {p,)Y u V{k)D ++ {k)G ++ {p,-q)G + -{p A ) 

The contribution of this diagram corresponds to a vertex correction in the ph channel by 
an effective interaction, which incorporates an intermediate particle-hole collective mode 
excitation. The contribution of the Fig. 2c reads 



i=l L 



d 4 Pi 
(2k)' 



dk 

(M 



(27r) 8 (5(g + p 4 -k- p 3 )5(k + p 2 - pi) 



Tr[r fl G-(q + p 4 )V(k)D-(k)G- + (p 3 )V(k-q)G + -(p 4 ) 
Tr [v{k)G- + (p 1 )Y u G ++ (p l - q)V(k - q)D ++ (k - q)G + '(p 2 ) 



(40) 



The latter diagram may be interpreted as a particle-hole fluctuation. The diagrams a-c are 
evaluated in the Appendix B. There we show that (i) the vector current contributions from 
diagrams a and b mutually cancel; (ii) the diagram c does not contribute because the axial- 
vector contribution involves traces over odd number of a-matrices and the vector-current 
contribution is cancelled by an equal and of opposite sign contribution from the diagram 
generated from c by flipping one of the loops upside-down; (iii) all contributions due to the 
Fermi-liquid interaction cancel after summing the diagrams a and b. For the contraction of 
the trace of the neutrino current with the polarization function we find = 1 ... 3) 



C d h(q, <?i,<?2) = iTr(Av) 



(?) + n i : + ' 6 (g) 



uiuj 2 - 




d A k 



G—(uYD—(k) 



(27T)*- 

G- + ( Pl )G + -(p 2 )G- + (p 3 )G + -( P4 ) 



(2ir) 4 5(q + p 4 -k-p 3 ) (2n) 4 S(k + p 2 - Pl ) . 



(41) 



This result is valid in the soft-neutrino and non-relativistic baryon limits. The second term 
on the r.h.s. in the square bracket can be dropped, as it does not contribute after the phase 
space integrations. Note that the total number of diagrams of the type a-c is four, if one 
allows for all possible relabelling of incoming and outgoing (identical) baryons; this forfactor 
is equal to the symmetry factor by which the total rate must be reduced. We do not include 
these factors explicitly. 



(39) 



D. Exchange contribution to the polarization function 

The exchange diagrams are generated from the direct ones by means of interchanging the 
outgoing propagators in a strong vertex. There is a complete set of diagrams analogous to a 
and b with exchanged labelling of the hole propagators. These contribute to the contraction 
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C cx (g) = lQg 2 A G 2 (J^j u^G" {uf J dkk 4 D~{kf 



r 4 

/n 



d 4 Pi 
(2tt)< 



G- + ( Pl )G + -(p 2 )G- + (p 3 )G + -( P4 ) 
{2n) 4 5(q + k + p 2 - p 3 )5{k + p 4 - Pi), 



(42) 



in the soft neutrino approximation. The skeleton diagrams which correspond to the inter- 
ference between the direct and exchange contributions are shown Fig. 3a-d. There are eight 
diagrams of each type if one allows for all possible relabelling of the propagators. 






Fig. 3: The exchange Feynman diagrams for baryon-baryon interaction in the 2p-2h approximation. Con- 
ventions are the same as in Fig. 2 

The analytical expressions for, e.g., the diagrams a and c are 



tt — h , a,ex 



(?) = u 



1=1 



d 4 pi 



dk 

(M 



■dk'(2ir) 8 5(q + p 4 - k - p 3 )5(k' + p 2 - p 3 )5(k + p 2 - Pi) 



Tr[r M G-( ? + p 4 )V(k)D-(k)G- + (p 3 ) 

V(k')D+ + (k)G + -(p 2 )V(k')G- + (p 1 )V(k')G ++ (q + p±)r u G + -(p 4 ) 



ill 



— h , c,ex 



(?) = u 



d 4 pi 
(2tt)< 



dk 



■dk'(27r) 8 8(q + p 4 -k- p 3 )S(k' + p 2 - p 3 )5(k + p 2 - Pl ) 



Tr[r,G-(q + Pi )V(k)D-(k)G-+(p 3 ) 

V(k')D ++ (k)G + -(p 2 )V(k')G- + (p 1 )T u G ++ (q + p i )V(k')G + -(p^ 

and their computation is a complete analogue of that for the direct diagrams (Appendix B). 
The vector current contribution again cancels among the diagrams a and c and, similarly, b 
and d. The contribution from the interference between the direct and exchange diagrams to 
the contraction of neutrino and baryon currents is 

CintO?) = l6g 2 A G 2 (J^ u^G—itu) 2 J dk J dk'k 2 k' 2 D—(k)D—(k') 



\ d l i 

i=i 



(2tt)< 



G- + ( Pl )G + -(p 2 )G- + (p 3 )G + -(p,) 
(2%) A 5(q + p 4 -k- p 3 )5(k + p 2 - Pi)5(k' +p 4 - Pl ), 



(45) 
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where we dropped the terms which vanish in the phase space integrations. The phase space 
integrations in the exchange contribution is complicated, since the momentum integrations 
do not decouple into two separate loops. The disentanglement can be achieved by constrain- 
ing the momentum transfer in one of the pion propagators at the value \k'\ = 2pp, as the 
main contribution to the integral originates near this value of the momentum transfer. 



IV. QUASIPARTICLE WIDTH 

The purpose of this section is to specify the width of baryon propagators. To this end we 
carry out a full resummation in the particle-particle (pp) channel by solving the scattering 
T-matrix at finite temperatures. Our approach is based on the Brueckner theory with the 
continuous energy-momentum spectrum of baryons. The non-perturbative treatment of the 
pp channel is mandatory for including the effects of the short-range correlations due to the 
repulsive part of the nucleon-nucleon force. These correlations are then responsible for the 
width of quasiparticle propagators, 7, in our perturbation expansion in the particle-hole (ph) 
channel. The ph interactions are dominated by the weaker long-range part of the nucleon- 
nucleon interaction, which makes possible the perturbative treatment of this channel by a 
truncation at two loops. The contour ordered T-matrix in the configuration space is: 

,%2] x 3, X A) — Y_fyp](xi, X2', X3, X4) 

+iV ipp] (x 1 ,x 2 ;x 3 ,x i )G(x 7 ,x 5 ) G(x & , x 6 )T(x 5 , x 6 ; x 3 , x 4 ), (46) 

where Vw,i(xi, x 2) x 3, x i) — c"z V[ pp ](xi, x 2 ; x 3 , x 4 ), is the time-local baryon-baryon interac- 
tion in the particle-particle channel. Note that the time locality implies that the pp propa- 
gator product GG = Gi™,] should be considered as a single matrix. The components of the 
scattering amplitudes, needed for complete specification of the self-energies, can be chosen as 
the retarded/advanced ones; the remaining components are provided by the optical theorem. 
In the quasiclassical limit the retarded/advanced T-matrices obey the integral equation 

T R / A (p,p';P) = VUPiP')+i J ^VUP,PlG^(p'',P)T R / A (p'',p',P), (47) 

where we kept the leading order terms in the gradient expansion of the product G [pp] R/A T R ' A . 
Here the subscript [pp] indicates the particle-particle channel and p, P are the relative mo- 
mentum and total four-momentum respectively. The two-particle Green's function, appear- 
ing in the kernel of equation fl4"7]), is defined as 



GBffe..*) = (ft/2 +R ) G> (ft/2-*; 



(2*)< 

-G< (P 2 /2 +Pl )G< (A/2 - Pl ) } W g^fij , (48) 

where we dropped the irrelevant dependence of the quasiclassical functions on their center- 
of-mass space-time coordinates. If the particle-hole symmetry is kept in the kernel of the 
integral equation, the T-matrix diverges at the critical temperature of the superfluid phase 
transition. To be able to apply our computation to the low-temperature regime (and thereby 
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avoid the pairing instability in the T-matrix) we drop the hole-hole propagators. This is 
a common approximation in the Brueckner theory and is justified in terms of the Bethe- 
Goldstone hole-line expansion. We treat the intermediate state two-particle propagation in 
the quasiparticle limit. Using the angle averaging procedure for the pp propagator and after 
partial wave expansion, the thermodynamic retarded T-matrix is given by 

Tg?(p,i/,P,u) = Vg liW (p,j/) 



+ -E /^/ 2 ^ p] ^(P,/)(Gg p] (/,P^))T^(/,p',P^), (49) 

7T v , J 



where a collectively denotes the quantum numbers (S, J, M) in a particular partial wave, p 
and P are the magnitudes of the relative and total momentum respectively, V(p, p') is the 

'\ppV 



bare nuclear interaction. Here (GfL-,) is the angle averaged two-particle propagator 



\G\pp\ [P, P, w)) - J — u _ e{P/2+p) _ £{P/2 _ p)+l5 . (5°) 

with e(p) = e p + 3?eS(£:p,p), i.e. the intermediate state propagation is treated in the 
quasiparticle approximation. The retarded self-energy is given by 

E R ( P ,u) = ±£(2J+1) Jdp'p n T^(j>,p>-p,p'-u + e(p'))f N {e(j>')\ (51) 

la 

which also defines its real and imaginary parts. The coupled equations ( [49"D and ([H]) are 
subject to normalization to the total density at a given temperature. 



V. PHASE SPACE INTEGRATIONS 



Let us turn to the task of evaluating the phase space integrals in the expressions for 
the current contractions. We substitute the Kadanoff-Baym ansatz in eq. fl4l|) and use 
the identity iW(£i)/iv( — £2) = 9(^1 — £2) [/at (£2) — In^i)], which is exact in the equilibrium 
limit. We then find that the contributions from each loop decouple, i.e., 

C div (q) = 16 9 2 A G 2 ( £l) u x u 2 G-(u>Y [^Lk*D--(k) 2 g(tu k )g(u J -u Jk )L(k)L(q-k), (52) 
\m^J J (2tt) 4 

where uj k = k and the elementary loop is defined as 

L(k) = ^ a( Pl )a(p 2 )[f N (e 2 ) - f^W^k+p, - Pl ). (53) 

The exchange contribution C ex leads to additional factor of two. The interference contribu- 



tion decouples only under certain constrains. The single loop, eq. (|53|) , can be evaluated to 



arbitrary order in the spectral width in general [2C]. We shall restrict to the small quasi 



particle damping limit and use the the expansion with respect to the width of the spectral 



function given by eq. (|30D . 
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1. Leading order 



The lowest order approximation corresponds to the quasiparticle (i.e. zero-width) limit. 
The contribution from a single loop vanishes in the time-like region of the phase space where 
uo k > \k\. This result is found only if the relativistic kinematics is applied; non-relativistic 
kinematics leads to spurious terms oc m/q. In the space-like region of the phase space the 
result is finite. We carry out the energy integrations keeping only the leading order term. 
Removing one of the trivial momentum delta functions we find 



L (k) 



d 3 p 



LM £ p) - fN(£ P +k)](2n)S(uj k + 



The integrations can be carried out exactly 
d 3 p 



L (k) 



[./V(£p) - f N (s p+k )]{2T[)5(uj k + e p -e 



11) 



* 2 



where m* is the effective mass of a quasiparticle and 

1 + exp [-p (e-(k) - /J,)) 



C(u k , k), 



C(u k , k) = In 



1 + exp [-p (e+(k)-fi)] 



(54) 



(55) 



(56) 



with e±(k) = {oJ k + e|)/4£fc ±uj k /2. Note that the quasiparticle loop (|55|) is zero in the time 
like region (u k > \k\), which sets a natural cut-off in the phase space integrations below. 



2. Next-to-leading order 

The next-to-leading order contribution (which is linear in 7) is 

/d^p f did 1 

—-^(27r)8(e + u k -e p+k ) J — 7 (j/) — 

x {8(e -e p )-S(e- a/)} [f N (e) - f N (e + uj k )) 



V 



(57) 



where we summed the two term arising from the product of the leading and next-to-leading 
order contribution to G < {p). The angular integral can be carried out analytically to the 
accuracy 0{ 7 2 ). One finds 



Am* 2 f ds 

Li(u k , k) = — J j—j [f N (e P ) - fN{s P + u)\ {Z(e p , k) - T[e v , k,u k )}, 

where the first term in the curly brackets is due to the wave-function renormalization 

V (u k - e q f 



Z{£p, fe) — &(Sp ^min) 



duj 7 {uo) 



1 no) ^min ■ 

{u - e p ) 2 

The second terms is the off-pole contribution and is given by 



4e„ 



J-{uj k , k, e p ) = arctan 



e/c — — H + 2^e p e k 
r/(e P + uj k )/2 



arctan 



£k — u k — [i — 2 y/ /e p e k 
7 (£ p + u k )/2 



(58) 



(59) 



(60) 
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The current contraction, which so far includes contributions to all orders in 7, now can 
be decomposed in the leading and next-to-leading order terms with respect to 7, employing 
the corresponding decomposition for the loops. E.g. for the direct contribution one finds 

g(u k ) g(u - u k ) [L (k) L (q - k) + Li(fc) L (q - k) + L (k) L x {q - k)\ . (61) 
The exchange and interference terms can be decomposed in a similar manner. 



3. Neutrino emissivity 

After the preparatory work above, the computation of the neutrino emissivity is straight- 
forward. We first relate the current contraction to our original expression for the neutrino 
emissivity by using the relation — 2g B (q )^smH^ l ,(q) = zll^ (g). Expression pi]) takes the 
form: 

*- = S / / wfei) / W?W S ^ + * - *> + C(?) < (62) 

where C(q) is the sum of the direct, exchange and interference contributions. Let us first 
compute the contribution from the direct term by substituting eq. ( 4lD for the current 
contraction. We carry out the integrations over the neutrino phase space and the summation 
over the three neutrino flavors to find: 
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-9 2 A G F 2 (—) [°° dujuj & G-~(ujf ! dkk 6 D—(kf 
\m n J Jo J 



uv 5(2tt) 7; 

J du k g{u k ) g(uj - u k ) [L (k) L (q - k) + L x {k) L (q - k) + L (k) Lt(q - k)} , (63) 

where we used d A k = Airk 2 dkdu k . Normalizing the energy scales by the temperature and 
the momenta by 2pp we obtain 

= 5(Sp G ^ (£) ' »"* = 5 ' 5 X 1019 /3T » <"* Cm " < M) 

where Tg is the temperature in units of 10 9 K, ^3 is the integral / in units 10 3 defined asQ 

roc roc roc ( 

1= dyy 6 G-(y) 2 Q(y) dxx A D (xf / dz g{z) g(y - z) {C(z,x)C(y - z, x) 

JO Jo J-00 [ 

2 2 1 

+ -zC(y - z, x) [F{z, x) - Z(z, x)} + ~{y - z)C{z, x) [F{y - z, x) - Z(y - z, x)} \. (65) 

7T 7T 



It is understood that the functions of new variables are relabelled. 
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The explicit dependence of eq. (|64"D on the temperature and density is the generic one • 
Additional dependence on these parameters is contained in the integral (|65D . Note that to 
avoid spurious contributions from the quasiparticle part, (the first term in curly brackets 
in (|65D ), one should restrict the z integration to the space-like region. For the numerical 
evaluation of the neutrino emissivity we use, following ref. H, the free space pion propagator: 



D (k) = [k 2 + m 2] 

The free-space approximation should be valid in the vicinity of the nuclear saturation density. 
The softening of the one-pion exchange (a precursor of the pion-condensation) increases the 
neutrino emissivity by large factors |J. We do not attempt to accommodate this effect as 
our main interest here is the role of the finite width of quasiparticles. The Pauli blocking 
factor 

Q(y) = 30 f'dww 2 (l-w)*[l-f v (wy)][l-M(l-w)y)], (66) 
Jo 

accounts for the occupation of neutrino and anti-neutrino final states. In the dilute (anti- 
neutrino limit /3fjL Uf <C 1 (where \i Vj is the chemical potential of neutrinos of flavor /) 

G(y) = l. 

In the low-temperature limit C(z) = z and the ^-integration decouples from the re- 
integration. On imposing 7(0;) — > (quasiparticle limit) one finds that T = and G (u) = 
oj~ 2 . Then the z integration can be carried analytically upon dropping the wave-function 
renormalization contribution: 

3 v (v 2 + 4n 2 ) 

dz g(z) g(y - z) z (y - z) = g ^ _ 1) ' (67) 

After these manipulations eq. (|64]) reduces to Friman and Maxwell's result [ref. [|J, eq. 
(47)]. The numerical coefficient in eq. (51), however, is by a factor 3 larger, since Friman 



and Maxwell do not carry out the summation over the three neutrino flavors at that stage. 

The contribution from the exchange current contraction, eq. (|42|) , leads to a factor of 2 
in the integral (|65|) . The contribution of the interference term, in the approximation where 
one of the momentum transfers is fixed at the characteristic value 2pp, is 

pOO poo poo ( 

hnt= dyy 6 G"(y) 2 Q(y) dxx 2 D" (x)D- (I) dz g{z) g(y - z)\ C(z,x)C(y - z,l) 

JO JO J~oo [ 

+ -zC(y - z, x) [F(z, 1) - Z(z, 1)] + -{y- z)C(z, 1) [F(y - z, x) - Z{y - z, x)} ) . (68) 

7T 7T 



VI. RESULTS 

The numerical calculations were carried out for pure neutron matter using the Paris NN 
interaction keeping J < 4 partial waves. Fig. 4 displays the real part of the on-shell self- 
energy and the half width of the spectral function as a function of the particle momentum for 
several values of the temperature at the saturation density n s = 0.17 fm~ 3 . The width of the 
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quasiparticle propagators can be parametrized in terms of the reciprocal of the quasiparticle 
life time in the Fermi-liquid theory (damping of the zero sound): 

(69) 

where a is a density dependent phenomenological parameter. The parabolic dependence of 
the width on the frequency is justified for temperatures below 30 MeV in the range of the 
densities n s < n < 2n s . The quadratic dependence of 7 on the temperature breaks down at 
slightly lower temperatures. The value of the parameter a weakly depends on the density 
and is approximately 0.2 MeV -1 . 



7 



aT 2 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 
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Fig. 4: The real part of the on-shcll self-energy and the half-width as a function of particle momentum at 
the saturation density n s = 0.17 fm -3 for different temperatures; the zero temperature Fermi momentum is 
1.7 fm - 1 . 
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The emergent neutrino spectrum can be caracterized by their spectral function 

rl poo ( 

S(y) = G-(y) 2 Q(y) / dxx A D"{x) 2 / dz g(z) g(y - z) \C(z, x)C(y - z, x) 

JO J-oo [ 

2 2 1 

+ -zC(y - z, x) [T(z, x) - Z(z, x)) + ~{y - z)£(z, x) \T(y - z, x) - Z(y - z, x)\ > , (70) 

TT 7T J 

which is depicted in Fig. 5. The values of the integral are shown as a function of neutrino 
frequency at T = 20 MeV and the saturation density n s = 0.17 fm -3 in the limit of vanishing 
width (dashed line), including the leading order contribution in the width (dashed- dotted 
line) and full non-perturbative result (solid line). The energy carried by neutrinos is of 
order of uj ~ 5T in all three cases, as the peak in the spectral function is independent of the 
approximation to the width of the propagators. The integral I3 is show in Fig. 6. The finite 
width of propagators leads to a suppression of the bremsstrahlung rate as a result of the LPM 
effect. Keeping the full non-perturbative expression for the causal propagators enhances the 
value of the integral, as the higher order terms contribute additively to the leading order 
result. The LPM effect sets in roughly when u ~ 7. As neutrinos are produced thermally, 
the onset temperature of the LPM effect is of the order of 7. Equation ( |5§D shows that the 
value of the parameter a controls the onset temperature which turns out of the order of 5 
MeV in agreement with the previous results of refs. |)7HlO| and our numerical computation 
(see Fig. 6). 




Fig. 5: The neutrino spectral function ( |70| ) at the temperature T = 20 MeV and density n s — 0.16 fm~ 3 . 
The dashed curve is the zero width limit, the dashed-dotted curve includes only the leading order in 7 
contribution from the causal propagator, the solid curve is the full non-perturbative result. 
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Fig. 6: The integral ( p5[ ) (including the exchange terms) as a function of temperature at the density 
n s = 0.16 fm~ 3 . The dashed curve is the zero width limit, the dashed-dotted curve includes only the leading- 
order in 7 contribution from the causal propagator, the solid curve is the full non-perturbative result. 

VII. CONCLUSIONS 

In this work we formulated a transport theory for neutrinos in the framework of real- 
time Green's functions formalism, with particular attention to the collision integrals for 
the neutrino-pair bremsstrahlung. The main focus was a first principle calculation of the 
bremsstrahlung emissivity including the width of propagators. This allows to answer ques- 
tions, not covered by the semi-phenomenological theory, such as the magnitude of the con- 
tribution of higher order terms in the expansion with respect to the quasiparticle width or 
the cancellation of the vector current contribution at all orders in the quasiparticle width. 
Even though the expression for the emissivity, which follows from our quasiclassical trans- 
port equation, is the same as the one found in the linear response theory, it is valid under 
conditions arbitrary far from equilibrium. This is particularly important in the regime where 
the neutrinos decouple from matter and their distribution function strongly deviates from 
the Fermi-Dirac form. 

The central quantity of the theory is the particle-hole polarization tensor in the ph chan- 
nel truncated at two loops. The pp channel is treated non-perturbatively within the finite 
temperature Brueckner theory. We find that the only contribution to the bremsstrahlung 
rate comes from the contraction of the tensor force with the axial vector current to all or- 
ders in the quasiparticle width. Other contributions, which arise from the contraction of the 
Fermi-liquid type interaction with the axial vector current and the contraction of the net 
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strong interaction with the vector current, cancel when we take the sum of the diagrams cor- 
responding to vertex corrections and propagator renormalization in the ph channel. Thereby 
the vector current conservation is established at all orders in the quasiparticle width. These 
cancellations are independent of the approximations to the propagators and are effective 
both in the quasiparticle limit and beyond. The three ingredients crucial to the cancella- 
tions are: (i) the anti-commutation of the tensor force with the axial vector current, (ii) the 
odd parity of the causal propagator under the exchange of its energy argument, (iii) the soft 
neutrino and non-relativistic baryon approximations. 

Our numerical evaluation of the neutrino emissivity of hot neutron matter, carried out at 
two loops, shows that the LPM-type suppression sets in at temperatures T > 7, in agreement 
with the previous work limited to the first order terms in the quasiparticle width (see ref. 
H and references therein). The higher order terms enhance the magnitude of the neutrino 
emissivity compared to the leading order result. The non-perturbative result, however, is 
still suppressed as compared to the quasiparticle limit. 

Our formalism can be extended in various ways. One obvious extension is allowing for 
two different chemical potentials of scattering baryons. This will include the Urea pro- 
cess (the /3-decay in the second order in the virial expansion) and the effects of the Pauli 
spin-paramagnetism, which become important in strong magnetic fields. The formalism 
can be adapted, with minor changes, for a computation of the space-like analogous of the 
bremsstrahlung and, in particular, the neutrino opacities of the supernova matter. The per- 
turbative scheme, employed here, itself requires further improvements in several direction, 
numerically the most important one being the renormalization of the one-boson exchange 
interaction in the ph channel. 
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The six Green's functions of the non-equilibrium theory are not independent. For com- 
pleteness we summarize here the linear relations among them, which can be easily verified 
from their definitions. The four components of the matrix Green's function are related to 
each other by the relations 
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APPENDIX A: REAL-TIME GREEN'S FUNCTIONS 



S (xi,x 2 ) = 6(t 1 - t 2 )S + ~(x l ,x 2 ) + 0(t 2 - ti)S' + (x l ,x 2 ), 
S ++ (xi,x 2 ) = 9(t 2 - t 1 )S + ^(xi,x 2 ) + 9(ti - t 2 )S~ + (xi,x 2 ) ) 
S (xi,x 2 ) + S ++ (xx,x 2 ) = S~ + (xx,x 2 ) + S + ~(x 1} x 2 ). 



(Al) 
(A2) 
(A3) 



Following Hermitian conjugation relations hold: 
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S (xi,x 2 ) = -S ++ *(x 2 ,xi), 
S~ + (x l} x 2 ) = -S~ + *(x 2 ,xi), 
S (xi,x 2 ) = -S + ~*(x 2 ,x l ). 



(A4) 
(A5) 
(A6) 



The retarded and advanced Green's functions are related to the components of the matrix 
Green's function via the relations 



S R (x 1 ,x 2 ) — 6(ti — t 2 ) S + (xi,x 2 )-S + (xi,x 2 ) 

= S~~(x lt x 2 ) - S~ + {x u x 2 ) = S + ~(x 1 , x 2 ) - S ++ (x 1 , x 2 ), 

S A (x 1 ,x 2 ) = -6(t 2 - ti) S + ~(xi,x 2 ) - S~ + (x, y) 

= S (xi,x 2 ) - S + ~(xx,x 2 ) = S~ + (xi,x 2 ) - S ++ (xi,x 2 ). 

They are Hermitian conjugates, i.e. 

S A (x ll x 2 ) = S R *(xi,x 2 ). 

In addition, we note that in the momentum representation they satisfy the equations 



S-{u,p) 



, S A (u,p)= S R (u,p 



(A7) 
(A8) 

(A9) 
(A10) 



The relations above are valid for the baryon and pion propagators in general, and we do not 
repeat them here. 

Similar relations hold among the self-energies. These can be identified by performing 
a unitary orthogonal transformation affected by the matrix R — (1 + ia y )/2 by means of 
formula 5" = R'^^SR. The form of the original Dyson equation in the matrix form (|3|) is 
invariant against the transformation, 



SL(xi,x 2 ) =£u)(xi,x 2 ) + S! J0 (x 1 ,x 3 )Qf_(x 3 ,x 2 )S!_(x 2 ,x 1 ) 
where the primed quantities have the "triangular" form 



q, 



12 



^ ^12 1 rV ^10 
gR gK } ! ^12 



'12 "12 



L 12 " 12 

n A o 



where 



S K (xx,x 2 ) = S c (x 1 ,x 2 ) + S a (x 1 ,x 2 ) = S > (x 1 ,x 2 ) + S < (x 1} x 2 ), 



(All) 



(A12) 



(A13) 



Q R {x u x 2 ) = VL c (x u x 2 ) +tt < (x 1 ,x 2 ), Q A {x 1 ,x 2 ) = Q c {x u x 2 ) + Q > (x 1 ,x 2 ), (A14) 



(xi,x 2 ) = Q c (x 1 ,x 2 ) + Vt a (xi,x 2 ) = -fP(xi,x 2 ) - fi < (xi,x 2 ). 



(A15) 
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APPENDIX B: DETAILS OF THE COMPUTATION OF THE POLARIZATION 

FUNCTION 



As an example we compute here the direct contribution to the polarization function, 
represented by the diagrams a and b in Fig. 2. The cancellation among the various contri- 
butions from these diagrams does not depend on the details of the structure of the baryon 
propagators (quasiparticle or dressed), but solely on the odd parity of the causal Green's 
function with respect to a change of the sign of the energy argument in the soft neutrino 
approximation. 

In the first step we substitute the vertices. As the contribution of the Landau Fermi- 
liquid part of the interaction will cancel out, to save space, we shall drop its contribution 
from the outset. For the diagrams a and b (excluding the factors for the topologically 
equivalent diagrams) we find 



*n- + '°(<z) 



' G 
2y/2 



Tr 



Tr 



Pi ft 



i=l 



d 4 pi 



dk 



(2tt) 4 J (2k) 4 
(<^0 - QA^m) G~~(q + p 4 ) (cr ■ fe) D~~(k) 

G- + (p 3 ) (cr . fe) D ++ (k)G ++ (q + p 4 ) (6 u0 - g A S vj ai) G + -( Vi ) 
o- • fe) G- + ( Pl ) (cr ■ fe) G + -(p 2 )} (2n) 8 6(q + p 4 - k - p 3 )5(k + p 2 - Pl ), (Bl) 



ill 



fJ,U 



{2V2) (X) /s 



d 4 pi 

WT 4 



dk 



Tr 



Tr 



(V - gA^iUi) G (q + p^ (cr ■ k)D (k) 

G- + (ps) (<U - 9aK 3 o 3 ) (cr ■ fe) D++(k)G++(p 3 - q)G+-(p i ) 
"(cr ■ fe) G- + { Pl ) (cr • fe) G + -(p 2 )} (2n) 8 6(q + p A -k- p 3 )5(k + p 2 - Pl ). (B2) 



Next we apply the approximation (j33|) to the causal and acausal Green's functions and fix 
their momenta at the corresponding Fermi momentum. Combining the diagrams a and b, 
we find 



ill-+< a (q)+ill-+< b (q) 



G 



m-n 



n 

i=l 



d 4 pi 
(2k)' 



dk 



G-(Lo) 2 D-(k) 2 G- + ( Pl )G + -(p 2 )G- + (p 3 )G + -(p 4 ) 

Tr{(5 M o - 9A^i(Ti) (cr ■ fe) (cr • fe) (5 v0 - g A 5 vj aj) 

- - 9A^i<Ji) (cr • fe) (5 v0 - gA^ujO-j) (cr • fe)} 

Tr [(cr ■ fe) (cr • fe)] (2n) s 5(q +p 4 -k- p 3 )5(k + p 2 - Pl ) 



(B3) 



where we used the conjugation relation ( |A10| ). The terms under the trace oc 5 0ai , 5q v vanish. 
The n o component of the polarization is hence zero and the vector current is conserved. 
The remainder simplifies to 
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m-n 



n 



(27T)' 



[G-(^) 2 J D-(fc) 2 G- + (p 1 )G + -(p 2 )G- + (j93)G + -(p 4 ) 



Tr 



(Jj (<T • fe) (<7 • fe) (Tj — CTj (<7 ■ fe) (7j (<T • fe) Tr [(<7 ■ fe) (<T ■ fe)] 



(27r) 8 (5(g +p 4 - * - Ps)S{k + p 2 - pi). 
The computation of the trace using the cr-algebra gives 
Tr [(tr • fe) (cr • fe)] Tr (cr • fe) (cr • fe) cr, — cr^ (cr • fe) (cr • fe) 



(B4) 



8k 2 - kikj) . (B5) 



The contraction of the polarization tensor with the trace of neutrino currents, given by 



TrAy = 8 [quq 2 j + gij<?2i + (^1^2 + <?i • c/2) ^ + e injm qlq^} 



leads to 



8/c 2 TrAy - fc^) = 128A; 4 



u 1 u 2 



(gi • fe)(gi • fc) 

fe 2 



(B6) 



(B7) 



Combining eqs. ( |B4| ) and ( |B7[ ) we recover eq. (ffl]). 

Let us turn to the fluctuation diagram in Fig. lc. From the original diagram one can 
generate three additional ones by turning each of the loops upside-down. Let us combine the 
diagram in Fig. lc with its counterpart, say d, which results from c by turning the upper 
loop upside-down. The analytical expression corresponding to their sum is 



n 



i=i 



d A Pi 

(27T)' 



dk 

(M 



G-{u) 2 D-{k) 2 G- + {p l )G + -{p 2 )G~ + {pz)G + -{p,) 



Tr 



(^0 - QA^m) (cr ■ fe) (cr ■ fe) Tr [(cr ■ fe) (5 u0 - QA^vj^j) (cr ■ fe)] 



Tr 



k) (S^o ~ 9A^i^i) (cr • k) Tr [(cr • fe) (S u0 - g A 5 uj aj) (cr ■ k)] 



(2n) 8 5(q + p 4 -k- p 3 )5(k + p 2 - Pl ), 



(B8) 



where we dropped q compared with k in the strong interaction vertex. The contribution due 
to the axial-vector current vanishes because the traces are over odd number of cr matrices; 
the contribution due the vector current cancels as these are identical for diagrams c and d 
and are of opposite sign. 
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